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METHOD AND APPARATUS FOR SEMI-BLIND COMMUNICATION CHANNEL 

ESTIMATION 

Field of the Invention 

5 

The present invention relates to multi-carrier communication systems, and more 
particularly to channel estimation in multi-carrier communication systems that employ 
orthogonal frequency division multiplexing (OFDM). 

10 Background of the Invention 

OFDM is a multi-carrier transmission technique, which divides available 
frequency spectrum of a communication channel into many carriers, often referred to as 
sub-carriers; and adjacent sub-earners are orthogonally phased to each other. Each of 

15 the sub-carriers is then modulated by a low rate data stream. As the sub-carriers are 
packed more closely than, for example, in frequency division multiplexing (FDMA), 
OFDM allows the frequency spectrum to be used more efficiently. In addition, OFDM 
does not require complex time switching, as in time division multiplexing (TDMA), and 
therefore does not suffer the overhead associated with time switching methods. 

20 FIG, 1 shows an OFDM system for transmitting and receiving information. A 

serial stream of data symbols is provided via input 101 to a serial to parallel converter 
102 that converts the single data stream to several parallel data streams. An inverse 
discrete Fourier transform (IDFT) module 104 processes the parallel data streams and 
produces a corresponding number of orthogonal modulated sub-carriers which are 

25 provided to a parallel to serial converter 106. In response, the parallel to serial converter 
106 provides a serial data signal to a cyclic prefix adder 108, and the cyclic prefix adder 
108 produces a transmit data signal at output 1 10, and the transmit data signal is 
transmitted on a communication channel. 

With further reference to FIG, 1, a corresponding received data signal on the 

30 communication channel is provided to a cyclic prefix remover 114 via input 1 12. The 
cyclic prefix remover 114 removes the cyclic prefix from the received data signal and 
outputs a single stream of data to a serial to parallel converter 116. Resultant signals 
from the outputs of the serial to parallel converter 1 16 are provided to a discrete Fourier 
transform (DFT) module 118, which provides a corresponding plurality of demodulated 

35 data streams to a parallel to serial converter 120. A serial data signal Is then provided by 
the parallel to serial converter 120 lo an equalizer 135 and to a channel estimator 130, 
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and the equalizer operates with the channel estimator 130 to determine the originally 
transmitted data from received data, and provide the received data via an output 135. 

The cyclic prefix is employed to address distortion in the communication channel. 
Adding the cyclic prefix comprises repeating the last few samples of each data symbol 

5 at its beginning, prior to its transmission. The length of the cyclic prefix should be 
chosen to be greater than or equal to the duration of the impulse response of the 
communication channel. This allows equalization of the channel distortion in the 
frequency domain by using a single tap scalar equalizer for each carrier, independently. 
However, in order to do this the response of the communication channel needs to be 

10 characterized. In practice only an estimate of the communication channel's 
characteristics is used, hence the need for a channel estimator. 

There are several methods of performing channel estimation, these include the 
following schemes; Pilot Symbol Assisted Modulation (PSAM), Blind Channel 
estimation, and a coded pilot method. Each of these is briefly described below. 

15 ps AM adds periodic transmissions of known symbols or pilots. Pilots comprise 

data that is known by both the transmitter and the receiver. Therefore, communicating 
pilot symbols allows the receiver to determine the difference between what was 
transmitted and what was received, and thus compensate for any variations in the 
received symbols that are caused by transmission between the transmitter and the 

20 receiver i.e. the communication channel. An estimate of the characteristics of the 
communication 'channel is required to provide such compensation across time and 
frequency domains of the communication channel. When the time and frequency 
characteristics of the communication channel are varying rapidly, as in mobile 
communication applications for example, channel estimation must be performed more 

25 frequently, hence the need for more pilots to be transmitted in order to maintain reliable 
communication. Thus, reducing the available bandwidth for data transmission. 

FIG. 2 shows a PSAM scheme graphically, where both pilot symbols 202 and data 
symbols 204 are shown in a three dimensional grid across time 206 and frequency 208 
axes, and where the vertical axis represents transmission power 210 of the data and pilot 

30 symbols. In PSAM, the pilot symbols 202 are inserted at intervals across time and 
frequency between the data symbols 204. Consequently, part of the signal energy and 
bandwidth of the communication channel is used for transmitting the pilot symbols 202. 
A received data signal in a PSAM scheme is passed through a 2-D Wiener filter, which 
essentially performs interpolation based on the statistics of the communication channel 

35 so as to estimate the characteristics of the communication channel between the pilot 
symbols 202, i.e. where the data symbols 204 are received. In this way, the data 
symbols 204 can be correctly recovered by taking the estimated time and frequency 
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characteristics into consideration to provide channel equalisation. One implementation 
of the PSAM scheme is in terrestrial transmission in digital video broadcasting (DVB- 
T). 

A PSAM scheme provides good channel estimation even when applied to time 

5 variant channels. However, when the normalized maximum Doppler spread is high, 
caused by fast changes in the communication channel characteristics, the frequency at 
which pilot symbols are requited increases in order to track such fast changes. This 
results in more bandwidth being required for pilot symbols, up to ten percent of the 
bandwidth of the communication channel, and leaving less of the bandwidth for data 

10 traffic- 
Blind channel estimation does not use pilots. Instead, the data symbols themselves 
are used to estimate the communication channel. Consequently, bandwidth of the 
communication channel is preserved. Several blind channel estimation schemes for 
OFDM are known, however, their tracking ability in a communication channel whose 

15 characteristics change or vary with time, Rayleigh fading time variant channels, for 
example, have not been as good as that of the PSAM scheme. 

The coded pilot method is described in US patent 5912876 by ETmimy where a 
main signal, comprising a quadrature amplitude modulated (QAM) version of a signal 
to be transmitted, and a pilot signal, are coded separately and transmitted as part of an 

20 OFDM signal. When the OFDM signal is received, the main signal portion is detected 
and an estimation of tht communication channel is determined from the detected coded 
pilot signal portion. Then the detected main signal and the estimation of the 
communication channel are used Lo estimate the signal that was transmitted. The coded 
pilot method is simple to implement, and the coding enhances the detection of the main 

25 and pilot signals, in the consequent channel estimation process. 

However, a transceiver using the coded pilot method is necessarily more 
complicated due to the coding in the transmitter, and detection of the codes in the 
receiver. In addition, a portion of the bandwidth of the communication channel needs to 
be allocated to support the transmission of the coded signals, thus reducing the usable 

30 portion of a predetermined ban d wi dth. 

Hence, there is a need for a channel estimation scheme that provides good 
performance in a communication channel having varying frequency and time 
characteristics, while preserving the usable bandwidth of the communication channel. 
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Brief Summary of the Invention 

The present invention seeks to provide a method and an apparatus for semi-blind 
communication channel estimation, which overcomes, or at least reduces the 
5 abovementioned problems of the prior art. 

Accordingly, in one aspect, the present invention provides a multi-carrier 
communication system comprising: 
a transmitter comprising: 

a data input for receiving an unknown data signal for transmission; 
10 a known data input for receiving a known data signal; 

a signal power ratio input for receiving a power ratio signal indicating a 
ratio for combining the unknown data signal and the known data signal; and 

a data combiner coupled to the data input, the known data input and the 
signal power ratio input, the data combiner for combining the unknown data signal and 
15 the known data signal in accordance with the power ratio signal to produce a composite 
output signal comprising discrete data signals* wherein each discrete data signal 
comprises at least a portion of the unknown data signal and at least a portion of the 
known data signal combined in accordance with the power ratio signal, the data 
combiner having an output adapted to provide the composite output signal to a multi- 
20 carrier transmitter, wherein the multi-carrier transmitter transmits a transmit signal on a 
communication channel, wherein the transmit signal includes the composite output 

signal; and 

a receiver comprising: 

a multi-carrier receiver for receiving a receive signal corresponding to the 

25 transmit signal on the communication channel, and the multi-carrier receiver having an 
output for providing a corresponding composite signal, wherein the corresponding 
composite signal comprises corresponding discrete data signals, and the corresponding 
composite signal being shaped by at least one signal shaping characteristic of the 
communication channel; 

30 a channel estimator having a known data input for receiving the known data 

signal, an input coupled to receive the corresponding composite signal, an input coupled 
to receive the power ratio signal, and an input for receiving at least one estimate of the 
unknown data signal, the channel estimator for estimating the at least one signal shaping 
characteristic of the communication channel from at least the corresponding composite 

35 signal, the at least the portion of the known data signal of at least some of the 

corresponding discrete data signals, the power ratio signal and the at least one estimate 
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of the unknown data signal, and the channel estimator having an output for providing at 
least one estimated communication channel characterising signal; and 

an equalizer coupled to receive the corresponding composite signal, the 
known data signal, the power ratio signal and the at least one estimated communication 
5 channel characterizing signal, the equalizer for configuring at least one of its signal 
shaping characteristics to compensate for the at least one signal shaping characteristic of 
the communication channel, the configured equalizer for shaping the corresponding 
composite signal accordingly, and the equalizer having an output for providing at least 
one subsequent estimate of the unknown data signal. t 
10 In another aspect the present invention provides a multi-carrier transmitting 

system comprising; 

a data input for receiving an unknown data signal for transmission; 
a known data input for receiving a known data signal; 

a signal power ratio input for receiving a power ratio signal indicating a ratio for 
15 combining the unknown data signal and the known data signal; and 

a data combiner coupled to the data input, the known data input and the signal 
power ratio input, the data combiner for combining the unknown data signal and the 
known data signal in accordance with the power ratio signal to produce a composite 
output signal comprising discrete data signals, wherein each discrete data signal 
20 comprises at least a portion of the unknown data signal and at least a portion of the 
known data signal combined in accordance with the power ratio signal, the data 
combiner having an output adapted to provide the composite output signal to a multi- 
carrier transmitter, wherein the multi-carrier transmitter transmits a transmit signal on a 
communication channel, wherein the transmit signal includes the composite output 
25 signal. 

In yet another aspect the present invention provides a multi-carrier receiving system 
comprising: 

a multi -carrier receiver for receiving a receive signal on a communication 
channel, where in the receive signal includes a composite signal, and the multi-carrier 

30 receiver having an output for providing the composite signal, wherein the composite 
signal comprises discrete data signals, and wherein each discrete data signal comprises 
at least a portion of an unknown data signal and at least a portion of a known data signal 
combined in accordance with a signal power ratio signal, the composite signal being 
shaped by at least one signal shaping characteristic of the communication channel; 

35 a channel estimator having a known data input for receiving the known data 

signal, an input coupled to receive the composite signal, an input coupled to receive the 
power ratio signal, and an input for receiving at least one estimate of the unknown data 
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signal, the channel estimator for estimating the at least one signal shaping characteristic 
of the communication channel from at least the composite signal, the at least the portion 
of the known data signal of at least some of the discrete data signals, the power ratio 
signal and the at least one estimate of the unknown data signal, and the channel 

5 estimator having an output for providing at least one estimated communication channel 
characterising signal; and 

an equalizer coupled to receive the composite signal* the known data signal, the 
power ratio signal and the at least one estimated communication channel characterizing 
signal, the equalizer for configuring at least one of its signal shaping characteristics to 

10 compensate for the at least one signal shaping characteristic of the communication 

channel, the configured equalizer for shaping the composite signal accordingly, and the 
equalizer having an output for providing at least one subsequent estimate of the 
unknown data signal. 

In still another aspect the present invention provides a method for determining 

15 received data in a multi-carrier communication system, wherein a received signal 

includes a composite signal received on a communication channel having transmission 
characteristics, wherein the composite signal comprises a plurality of discrete data 
signals spaced in time and frequency, and wherein each discrete .data signal comprises a 
data portion and a pilot portion, wherein the data portion comprises one of a 

20 predetermined group of symbols, the method comprising the steps of: 

a) defining a set of the plurality of discrete data signals; 

b) setting a predetermined number of iterations; 

c) receiving the composite signal on the communication channel; 

d) selecting one of the plurality of discrete data signals to be estimated; 

25 e) selecting a group of the plurality of discrete data signals, wherein each discrete data 
signal of the group of the plurality of discrete data signals is relevant to determining the 
one of the plurality of data signals; 

f) normalising the group of the plurality of discrete data signals using at least the pilot 
portion of the one of the plurality of discrete data signals; 
30 g) estimating the transmission characteristics of the communication channel using the 
normalised group of the plurality of discrete data signals in step (f); and 
h) estimating the data portion of the one of the plurality of discrete data signals using 
the communication channel having the transmission characteristics estimated in step(g) 
and at least the pilot portion of the one of the plurality of discrete data signals. 
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Brief Description of the Drawings 

An embodiment of the present invention will now be more fully described, by 
way of example, with reference to the drawings of which: 
5 FIG. 1 shows a prior art communication system; 

FIG. 2 shows a graphical representation of data and pilot symbols communicated 
on a communication channel of the communication system in FIG. 1; 

FIG. 3 shows a communication system in accordance with the present invention; 
FIG. 4 shows a graphical representation of data and pilot symbols communication 
10 on a communication channel of the communication system in FIG. 3; 

FIG. 5 shows a flow chart diagram of a process in a receiver of the 
communication system in FIG, 3; 

FIG. 6 shows a graphical representation of the APSB scheme in the 
communication system in FIG. 3; and 
15 FIG. 7 and 8 shows comparative simulation results of the performance of the 

communication system in FIG.3 and prior art communication systems. 

Detail Description of the Drawings 

20 The present invention combines the advantages of the PS AM and the blind 

schemes to produce a scheme where data and pilot symbols are combined prior to 
transmission, and separated when received at a receiver. In addition, separation of the 
pilot symbols at the receiver is accomplished by treating the data as noise and applying 
an iterative process to detect the data symbols. An embodiment of the present invention 

25 will now be described. 

In FIG. 3 a communication system 300 that communicates data using an added 
pilot semi-blind (APSB) scheme, in accordance with the present invention, comprises a 
combiner 305 having a data input 307 for receiving a data signal 306, a pilot input 309 
for receiving a pilot signal 308, and an input 311 for receiving a data-to-pilot power 

30 ratio (DPR) signal. The data signal 306 comprises a series of data symbols while the 
pilot signal comprises a series of pilot symbols. The DPR signal includes information 
on power levels of the data signal 306 and the pilot signal 308 that are to be combined. 

The combiner 305 includes a data conditioner 312 that receives the data signal 
306 and the DPR signal and provides a conditioned data signal having a power level as 

35 indicated by the DPR signal Similarly, the combiner 305 includes a pilot conditioner 
313 that receives the pilot signal 308 and the DPR signal, and provides a conditioned 
pilot signal having a power level as indicated by the DPR signal. 

i 
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The combiner 305 also includes an adder 314 for combining the conditioned data 
signal from the data signal conditioner 312 and the conditioned pilot signal from the 
pilot signal conditioner 313, by combining a series of conditioned data symbols and a 
series of conditioned pilot symbols, and providing data-pilot signals 310 via an output 

5 3 1 8 of the combiner 305. The data-pilot signals 310 comprises a series of discrete data 
signals, where each discrete data signal has a data portion and a pilot data portion, and 
where the power levels of the data portion and ihe pilot portion are in accordance with 
the DPR signal. The data portion represents at least one data symbol 

The data-pilot signals 310 are then processed, in sequence, by the serial to parallel 

10 converter 102, the IFFT 104, the parallel to serial converter 106, the cyclic prefix adder 
108, as described earlier* and a transmit data-pilot signal from output 310 is transmitted 
on a communication channel. 

A corresponding receive data-pilot signal 314 received from the communication 
channel at input 315, is processed sequentially by the cyclic prefix remover 1 14, the 

15 serial to parallel converter 116, the discrete Fourier transform 118 and the parallel to 
serial converter 120, as described earlier. An output data-pilot signal 316 from output 
317 of the parallel to serial converter 120 is then provided to an APSB equalizer 325 
and to an APSB channel estimator 320, The APSB channel estimator 320 aad the APSB 
channel equaliser 325, each include an input 319 for receiving the pilot signal 308, and 

20 each include an input 3 18 for receiving the DPR signal. The pilot signal 308 and the 
DPR signal may be stored in a memory (not shown). The APSB equaliser 325 operates 
with the APSB channel estimator 320 to determine the originally transmitted data 306 
from the data-pilot signal 316 using the pilot signal 308 and the DPR signal, and 
provides the originally transmitted data 306 via output 330 of the APSB equaliser 325. 

25 This is accomplished with an iterative process where a channel estimate h(n) i$ made 
by the APSB channel estimator 320, and a data- estimate is then made using the 
previously obtained channel estimate and provided via output 330. The output 330 is 
coupled to the APSB channel estimator 320 to provide the channel estimate for a 
subsequent iteration of channel estimation and data estimation. When a predetermined 

30 number of iterations have been performed, the originally transmitted data 306 is 
determined and provided via output 330 of the APSB equaliser 325. 

With additional reference to HG- 4, which shows the APSB scheme graphically, 
symbols of the data-pilot signal 317 are shown as discrete signals 405 in the three 
dimensional grid across time 206 and frequency 208 axes, and where the vertical axis 

35 represents transmission power 210. Each discrete signal 405 comprises two portions, a 
data portion 407 and a pilot portion 409. The proportions are defined by the DPR signal, 
which is, as defined earlier, the ratio of the data signal power to the pilot signal power. 
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There is an optimum DPR where the performance of the communication system 300 
would be best. 

In accordance with the present invention as described, the pilot is transmitted 
along with the data simultaneously and in the same frequency band, therefore 

5 transmission of the pilot advantageously does not consume bandwidth of the 
communication channel. 

The data-pilot signal 314 received from the communication channel contains 
information about the communication channel, and the pilot portion of the data-pilot 
signal is known. Then, by treating the data portion as noise, an initial estimation of the 

10 communication channel can be determined, and a data estimate obtained, with the initial 
estimation of the communication channel, using a least squares approach. With partial 
knowledge of the characteristics of the communication channel via the channel 
estimates and the data estimates, further iterations of channel and data estimation are 
performed, thereby improving the accuracy of the estimation and allowing an accurate 

15 estimate of the transmitted data to be made after a predetermined number of iterations. 

With additional reference to FIG. 5, a process 500 for estimating the transmitted 
data starts 505 with selecting 510 a set of cells S D from the received set of data-pilot 
cells S. Next the predetermined number of iterations is set 515, and a first cell, for 
example cell 430 at (nj) is selected 535. Subsequently, cells relevant to the estimation 

20 of the cell 430 are selected 540. Here, the selected cells are the eight cells adjacent to 
the cell 430 as typically the closest surrounding cells are selected. Using the time 206 n 
and frequency 20S / axes as references, the cell 430 is identified as cell (n,l) and the 
relevant selected cells as S p (n f l). All the relevant selected cells S F (n,l) are then 
normalised 542, for example by dividing the selected relevant cells by the pilot symbol 

25 portion 409 to allow a suitable channel estimate to be made. An initial channel estimate 
545 is then made using a 2-D Wiener filter, where the data portions 407 of the selected 
relevant cells are treated as noise, and the result is an estimate of the communication 
channel at (nj). Subsequently, a data estimate 550 is made using the least squares 
approach, and a decision 555 is made as to which symbol was transmitted. The decision 

30 step 555 may involve processes such as interleaving and coding before a decision is 
made. 

A determination 560 is then made as to whether all the cells in the set S have been 
determined. When not all the cells in the set S have been determined* then the next cell 
of the set S of cells is selected 570, and the process 500 returns to step 540 of selecting 
35 relevant cells to the newly selected cell. The looping back through step 570 continues 
until all the cells in the set 5 have been determined. When all the cells in the set S have 
been determined, the determination at step 560 is true, a counter (not shown) indicating 
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ihe number of iterations z is incremented 562, and a determination 565 is made as to 
whether the number of iteration i have reached the predetermined number of iterations 
set in step 515- When the number of iterations i has reached the predetermined number 
of iterations set in step 515, the process 500 ends 567. However, when the number of 
5 iterations / has not reached the predetermined number, the process 500 returns to step 
535 and repeats as described above. 

Returning now to FIG. 3 a functional description of .the communication system 
300 now follows. The rules for notation that are employed are: 

in the time domain, signals at index (n,l) are written* as h(n,l), vectors as h(n), and 
10 matrices as H(n); 

in the frequency domain a subscript F is added (e.g. hi(n,t) t htfn), Hj(n)); and 

sub-vectors/macrices are capped with a tilde e.g. h or H r and estimated entities 

with a hat e.g. h; 

superscripts \ H and * denote the operations matrix transpose, Hermitian and 

15 conjugation, respectively, and E denotes expectation; and 

elements of vectors/matrices are denoted as [h]/[H} iir 
The elementary sample period of an OFDM system model is T, the number of 
subcarriers is L, the number of cyclic prefix samples is D, and the maximum number of 
channel response samples is 

20 5+1 (<Z>) 

while the total number of samples in one OFDM symbol is 

and, the OFDM symbol period is defined as 

7>pr. 

25 The complex baseband representation of the communication channel 605, 

particularly a mobile wireless communication channel, impulse response at time t is 
described by 

ft(^)=Xy/(')5(T-T/o) (i) 
i ; 

30 where T/fr) and Y;(0 are the delay and complex amplitude.; of the I th path, respectively. 
The power delay profile of the channel is defined as 

p{X)^E[h{t^h\t^] j : (2) 

Assuming an exponentially decaying power delayjprofile with 
p(T)=Acxp(-T/T wl ) ! 
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15 



20 



25 



30 



where % rmS is a parameter of the channel known as the root mean squared delay spread, 
and A is a normalizing constant. The normalized T niu is given as X^IT. 

Due to the relative mobility of the transmitter and receiver, a maximum Doppler 
spread, / 0T will occur in the received signal. This is accounted for with a time-variant 
channel where a high/^, implies a fast varying channel. The normalized maximum 
Doppler spread is defined as/^7^ based on the assumption that the channel coefficients 
are time invariant over each OFDM symbol period T y 

A function which is useful for analysis in OFDM is the time-variant transfer 
function obtained from the Fourier transform of equation 
delay T, which produces the equation below. 



'1) above, with respect to 



h F {tfl = J h{tX) exp(-;27l/C) dl 



(3) 



Assuming a wide sense stationary uncorrelated scattering (WSSUS), and a 
Rayleigh fading channel with Jakes" spectrum, the autocorrelation of the channel is 
separable in time t, and frequency/, and can be written as 
r hfh ,(t,t';f,f)=l{h F (t,f)h;(?,f)\ 
= r [ (Af)r / (A/) 



where Ar = t- f; and Af=f-f'.'We have 

1- exp(-Z>7-(l/r rm , +J2J*f)) 



follows. 



(4A) 



r / (A/)= 



(l-exp(-Dr/T rmI ))(l + >27rT n 
1 

: 1 + ;2AA™ 



) 



(4B) 



assuming [{D)Hr„JT)) » I, and 

r f (Ar)=/o(27iFz)AO (4C) 

with 7 D (-) being the zeroth order Bessel function of the first kind. 
The discrete channel model of a communication channel of order B( < D) can be 
described as h(n,l) - h(t = nT s ,T = IT), assuming that r^is uniformly spaced at 
intervals of T and that p(*)= 0 for r ) = BT and represented as a vector at time 
nas h(n)=[h(n,0)..Mn 7 L-T)] T 

=[h(n,0)...h(n,B)0...0] T eq^* 1 ( 5 ) 
where C is a complex number. 
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The channel coefficients in the frequency domain are obtained in a similar way 
to equation (2), but in discrete time as 



h F (n)=VlTb (n-)=[h F (n,0)...h F (n 1 L- 1)] 



(6) 



5 where F is the L X L unitary discrete Fourier transform (DFT) matrix with 

[F]^ j = exp(-27OT//L)/-\/r, andF H is the corresponding inverse .DFT (IDFT) matrix. 
The auto-correlation of the discrete channel, r h h (n,Z)> can be obtained from the 
equations (4A) f (4B) and (4C) by replacing &f=I/(LT) ajnd At=nT s . 

The received OFDM signal may be considered a digital signal in a 2-D cell 

10 structure with indices (n,Z)eS- 

where S={{n,l):Q <n<N -1,0 <1<L-1} , and where N is the total number of cells in 
the time direction, n, and L is the number of sub-carriers as defined earlier. 

A model for a 2-D Wiener filter will now be described, as such a filter is used to 

i 

form estimates from sampled signals at the receiver. Estimates will be obtained at 
15 index (nJ)eS 0 from the sampled signals at indices (rf J )eS£ n, '\ where S a and S p n ' l) 
are subsets of S y and where S p =u n ^ S F n,l \ The definitions of S P ,S P andS P n,l) will be 
provided later j 

In a 2-D wide sense stationary (WSS) stochastic process h F (n 7 l) 

which contains information of the desired signal h F (n J)zydi$ corrupted by correlated 

i 1 

20 noise h F (nJ)v F (nJ) and additive white Gaussian noise (AWGN), u F (n,l), as follows: 



(7) 



25 



30 



35 



The notation h F (nJ) and u F (n,l) indicate that the channel frequency response 
is estimated from its noisy samples. It is assumed that v r (n, I) and u F (n,l) are 
white, and that /* F (n,/),v^(/i,Z)andu F (n,/)are mutually un-correlated, and zero mean 
stochastic processes. In the prior art PSAM scheme, the correlated noise component is 
absent, m contrast, for the APSB scheme of the present indention, as described, the 
correlated noise component is advantageously reduced through successive iterative 
measures. 

For the 2-D Wiener filter, the estimator for h F (n,tyis given as: 



(8) 



where w(n,l\riT ) is the weight applied on h F (riJ )to estimate h F {n,l\ 

Based on the assumption that a total of A N and A/j pilots are used in the n and I 
direction, respectively, for each estimation of h F (nJ) , a vector 6(n 7 /)eC A * Al x 1 is 
performed by stacking the elements from h F (ri ,i ), V(«' J )eS)> ' ; and a vector 
(n,f)eC A " AL X ^ formed by stacking the conjugate elements of w(n,/;« r /' ) . 
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Equation (8) can be re- written as follows. 

h F (nJ) = W*(n,0 h F (nJ) 



To minimize the difference h F (n,t) - h(n,l) in the mead square sense and obtain the 



optimum tap- weight vector, the orthogonality principle is 
5 resulting in Lhe equation below. 



(9) 



now applied to equation (9), 
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E [(M^0-V(^0)^(^^/ ,, )] = 0, V{n",r} e S^ l) (10 ) 
Substituting equation (8) into equation (10), the Wiener-Hopf equation, with optimum 
tap weights w 0 (nJ;n\O , is obtained as provided below, 

£ yv 0 Mn\nE[h P {n\nhp(n\n\ = 4h F (rt f[ 0/£(n'\n] (11) 
{n\f'} esp» 

Defining the cross-correlation of h F {nJ)andti F fyi,l) and the autocorrelation 
of h F (n,l), respectively as 



15 



r R {n-n\l-n=E[h F (nMi(r( 



20 



Letting. r /U (n,t)=E[ h F (n,l)h* F (n>l)]eC*** L * 1 be the cross- 
correlation vector and £ (h,/)=e[ h fC^O h?(nJ)]eC A "^* A * Ai be the 
autocorrelation matxiK h F (nJ) formed from the elements of equations (12) and (13), 
respectively. Therefore, the Wiener-Hopf equation can be re-writtcn in matrix notation 
25 as 



(12) 



(13) 



(14) 



Assuming that hpfafyvfafyanduinj) are mutually un-correlated and wide- 
sense stationary white stochastic processes, using the definition provided by equation 
30 (7), and letting An and AZ be the discrete time and frequency difference indices, 



respectively, equations (12) and (13) can be written as fol 

(An, AO = ^ F (An, AO + (An, AO x r Vf Vp (An, AO 



35 



-{ 



%h F (An.AO+cT^+cr^ 



lows. 



(15) 



forAn 7 A/ = 0 
otherwise 



(16) 
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15 



20 



25 



30 



35 



where tj ^(An,AZ) is the autocorrelation matrix h F 
R £ ^ and r /^( n >0 can be simplified as 



R 



(17) 

n r M >z)= r M >0 J (18) 

And, where R h F h F is the NL by ML autocorrelation matrix of h F% then R,^can be 
obtained from equation (4). If R /IfA/? o\and(j\ are known, the optimum weight is 

W 0 (n,0=[RA^+(^W+^5) i] 1 Th^M (19) 
The mean square error (MSE) can then be obtained as j 

a \/(n,l)eSn\ (20) 
(nj) is then! obtained by substituting 



\h F (n,l)-h F (n,l)\ 
The minimum MS# (MMSE) at index 
equation (9) into equation (20), and using the optimum weight factor in equation (19), 

(21) 



MMSE(n,l)=<J 2 k - w? fof)r Mf (n f 0 



For simplicity, tr /r =l for all subsequent equations. 
Now letting the transmitted data in the frequency domain be 
d F (n,l),\f(nJ)eS D > where n is the discrete time index |and Z the discrete frequency 

i 

index. In an OFDM system, n refers to the OFDM block index, while I refers to the sub- 
carrier index, and pilots are denoted as p F {n,l), V(n,Z)e Sjy. It is assumed that Pf(nJ) is 
deterministic and selected from a fixed set of alphabets, v/hile d F (n, I), the data is a 
zero mean stochastic process. The variance of p F (n,l) is 'denoted as s l o^ , and in 
practical implementations s l {(\ . j 

In accordance with the present invention, as described herein, pilot signals and 
data signals axe added together and co-exist at all time, /i,jand frequency points T /, in 
order to conserve bandwidth, that is, S D = S P -S- The signal after the data signal and the 
pilot signa] are combined is defined below. , 

x F (n, I) = d F (n,l)+p f (n.I), V(n, l)eS (22) 

and c x =(1h-£ )<v Other definitions follow. ; 

Data -to -Pilot Power Ratio (DPR)= a] cr^) = l/e 2 

Signal -to- Pilot Power Ratio (SPR)=a 1 z /(*j 2 cr^) = 1+ 1 / e 1 = 1 + DPR 
Assuming that the orthogonality of the OFDM system 600 is maintained, the 
signal after OFDM demodulation is 

y F (n,l)=h F (nJ)x F (nJ)+b F (nJ) 

=h F (nMd F (nJ)+ PF (nJ))+b F (n,l) (23) 
where h F (n,l)is an additive Gaussian noise (AWGN) stochastic process. 



10 



15 



20 



35 



15 

In accordance with the present invention, as described herein, channel-data 
estimation from the received signal y F {n,T) is performed injan iterative manner. For ease 
of description, the first iteration and subsequent iterations will be described separately. 

FIRST rrSKATTON j 

Channel estimator 320 is a 2-D Wiener filter, as characterized earlier, estimates the 
channel response. First, the received signal y F (n 7 l) in equajtion (23) is normalized in the 
first iteration, thus. 

h F (nJ)=y F (n,l)f p F (n,l) 



=h F (n,l)+h F (n t l)v Fi (n,l) + u n (nj) 



(24) 



where v n (nj) = dp(njy p F (n 7 l) is a data-dependent noise introduced due to the addition 
of the pilots to die data, and u n (n % I) = b r {nj) I p F (n,t) is AjiVGN. The subscript "1" 
indicates that the notation is specific to the first iteration. Similar use of subscript "i" 
will be employed for the iteration. 

Assuming that the piJots 7 p F (n 7 l) 7 are selected from) i 
symbols, we obtain 



and similarly 



2 

CF vl : 



2 



p F (nJ) 



_ <*d _ 



a set of constant modulus 



DPR 



e 2 cr?~SNR 



(25) 



(26) 



where, as defined earlier, SPR = (1 + e 2 ) I e 2 , and SNR = cr 2 



the optimum tap weights can 



Equation (24) is in a similar form as equation (7) and 1 
be obtained in a similar form as in equation (14) using time A N and frequency A^ 
25 sampled signals. In this case, the autocorrelation matrix s > i$ independent of (nj), 
however the cross correlation vector T^ hf is usually dependent of (n, I) . 

FIG. 6 shows a graphical representation in the time/frequency domain, and S P n,!) 
is selected such that (n f l) is in the center of S F n,l l The set of indices encompased by S P n ' l) 
will change, and may be visualized as a sliding window, as (nj) .changes. The relative 
30 index of (n.l) and ^"' /} however, will not change. Assuming wide sense stationary, 
R - h (nj) will also not change. Next, for example, the cross-correlation vector 
r h F l F ^ n ^ w hen estimating index 601 and 602 would be the same. When (n,0 is near to 
the boundary of S Ft such as when estimating index 603, the relative indices of (n t l) and 



V° would be different. Thus, r^OU) would be different in this case. 
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Now, the MMSE estimator for the first iteration follows. 



Thus, the first estimate of x F (nJ) is 

Then, in view of equation (22), the estimate for d F {n,l) is obtained as 
d n (nj)= x n (nJ)-p F (nJ) 



(27) 



(28) 



_ ft p (n, Z)fr*r (n, Q + p F (n, 7)) + ^(n, Q 



10 



h F[ (n t l) 



If the channel estimation is perfect, i.e. h F1 (nJ)= h F (n,l), 



Then, a decision device or slicer, as is known in the art, is 
an estimate 3 Fl (nJ) of d F (nJ). 



15 



SUBSEQUENT TTERATTQNS 
For the second iteration, the normalization is carried out using d F} (n,l) + p F (n,l). 
The normalization tries to remove the data-dependent : 
the first iteration as characterized in equation (24). The normalization for the second 
20 iteration is carried out as follows: 



-p F {nJ) (29) 



then equation (26) becomes 



(30) 



used on d n (n,l) to obtain 

i 



h Fn {n,l)--^ 



d Fi (n,0 + p F (nJ) 
h F fa l)(d F fa Q + p F fa, 1)) \ b F (n 7 -l) 



25 



d Fl {n,l)+p F (n,l) 



(31) 



Making the substitution d n (n,l) ~ d F (n 7 l) which is a good approximation when 
the probability of symbol error is small, equation (31 )becomes 

b F (n,l\ 



h F2 {n i l) = h F (nJ)+ 



d Fl (nJ) + p\r(nJ) 



30 where 



= h F (n,l)+u F2 (n,l) 
b P {n,l) 



(d n (nJ) + p F (nJ)) 



(32) 
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Equation (32) has the form of equation (7) with cr^=0 since v n (nj)=0. 

Then based on the assumption that £ 2 ({l 7 rf F (n 7 Z) is modulated by quadrature phase 
shift keying (QPSK), and p F (nJ)e{sG d (±l±j)/<>f2) , a QPSK constellation 
\d Fi (nJ)+p F {nJ)paxL have four possible values, each being equally likely to occur: 

^[(1+ j)+e(±l±jff2 . It should be noted that only he first quadrant of the QPSK 
signal constellation has been considered due to its 2-D symmetry. Hence, the required 



vanance reduces to 



(i-s 2 ) 2 (H-e 2 ) 



=(l+ £ 4 ) 



(l-2£ 2 +e 4 )(|l + £ 2 )tT^ 



*{l±2e 2 )SNIC x 
This approximation can be made when B 2 ((l, 



The 2-D Wiener filtering is applied accordingly k 



setting by setting 2 = 2 = 0 (i.e. there is no correlated noise component) and 



a l~ a li Then the channel estimator is 



Finally, the estimate d F2 (n,t)\$ obtained by following the steps as in equations 
(28) and (29). 

The subsequent iterations can be extended from the second iteration. However, 
the assumption d F2 (nJ)=d F (nJ) made in equation (32) will be more accurate in 
25 subsequent iterations and will therefore result in a better estimate for 
hf^n^.hf^nj), etc. As more iterations are performed, l 
converges to the actual channel response. 

In the prior art PS AM scheme, for every one pilot! < 



(33) 



equation (9) and (14) by 



0 



(34) 



cell, there are on average 
($n&l -1 ) data ce ^ s - The effective average signal power after pilot insertion is 



S L S N 



Lhe channel estimate 



where 7} 2 > 1 is lhe ratio of the power of the pilot jo that of the data for PSAM 



scheme. 



(35) 
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The signal-to-noise ratio (SNR) of the PSAM scheme can 



expressed as 



SMR=101og l0 (a2 /trl)=SNS! +SMR m {dB) 



SNR is the SNR of the system without the 

represents the additional SNR required to compensate for 
been allocated to the pilot. Therefore, it follows from 



(36) 

presence of pilots, while SNR )o5t 
the loss of power which has 
ions (35) and (36) that 



equati 



flV^=I01og 10 ((7j/o3)iB 



=101og 10 



7? 2 +(<5,5„-l) 



S L S N 



Since part of the bandwidth is used for pilot transit ssion, there is a reduction in 
the effective data rate, The percentage bandwidth loss is 

1 



on (36) for proper comparison 
additional power is incurred 



When power boosted pilots (i;e. 7))1) axe used, we 
actual SNR would have to be adjusted according to equati 
with other schemes. In the case when r/= 1 , however, no 
under the definition provided by equation (37). Regardlesjs of the value of 77, a 
bandwidth loss is still incuixed a$ seen from equation (38 

For the APSB scheme, in accordance with the preient invention, as described, 
the SNR is defined in a similar way as in equation (36), however with the APSB 
scheme 



thus be 



(37) 



(38) 



note that SNR,^ > 0 dB. The 



SNR lM =mog 1Q {(il/<J 2 d )dB 
r SPR } 



=101og 10 



{SPR-IJ 



dB 



In contrast to the prior art PSAM scheme, the 
greater than 0 dB since £) 0. On the other hand, although 
suffers from bandwidth loss, the APSB scheme advantag^ 

Both the prior ait PSAM and the APSB schemes 
symbols to be received before an optimum channel 
Depending on the selection of the indices of the sampled 
processing delays would be incurred in order for an optimum 



(39) 



AP$B scheme SNR j0SS is always 
the prior art PSAM scheme 
ously'has W jt>ffl = 0. 

t squire sufficient OFDM 
estimate can be obtained, 
signal, S/ 1,0 , different 
estimate to be made. 



i 
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When a non-causal 2-D Wiener filter is used, such that the index (nj) is at the center of 

the processing delay is then S s (AN- 1)/ 2 for the ririor art PSAM scheme. For 
the APSB scheme, the processing delay is the same with l' N = 1- 

The results of simulations will now be presented. For the simulations, the 
5 common system parameters were set as L = 32, D =8, B =7, t^ /T - 0.5, IT = 8 
and SNR = 20dB. It was assumed that the SNR is known and that a]- 1, and a QPSK 



signal was used In addition, 
of the 2-D indices of 
estimate the channel at 



signal constellation for the transmitted 1 signal and the pilot ; 
S P inJ> was selected such that the index is in the center t 
This results in the best MSE performance. It is possible to i 
10 indices 

n<(A N -l)/2orn>N-(ti N +l)[% and /<(A L -l)/2ori>L-(A L +l)/2 
in such a manner where a centered-inteipolation i$ carried out. For the other indices, an 
off-centered interpolation is required to estimate the channel, giving some degradation 
of the MSE, 

15 With reference to FIG. 7, consider first, the bit error rate (BER) performance at a 

maximum normalized Doppler spreader,, of 0.005. Setting A„J= & L = 15, the BER for 
the APSB scheme using different £(SPR) for the first to the fifth iterations, A high SPR 
means that low amount of power is used for pilot transmission while a small SPR means 
that low amount of power is used for data. Both extreme cases imply that h F (n 7 l) will 

20 be badly estimated. Hence, for each iteration, there exist an optimum point whereby the 
SPR wi)) minimize the BER. It is observed that as the number of iterations increases, 
the performance improves, but at a diminishing marginal unount. 

For the prior art PSAM scheme, as we increase T[, ' 

i 

improves. However, since SNR l0SE also] increases, less pows 

i 

25 a given SNR. Thus, the improvement of the BER brought ; 

estimation would be offset at some point as T) increases. Consequently, there is an 

optimum T| that minimizes the BER as well. 

i 

With both the prior art PSAM scheme and the APSB scheme, performance 
improves when larger filter taps are usjed i.e. larger A^ and A £ . For the PSAM scheme, 

30 the performance also improves for smaller A„ and A^ when closer sampling intervals 
are employed, while for the APSB scheme, more iterations (K) results in better 
performance, although the marginal performance gain would decrease for both cases. 
Hence, the system design is flexible since increased complexity of the PSAM scheme, 
and the APSB scheme would improve; performance. However, the APSB scheme, in 

35 accordance with the present inventionias described> advantageously does not suffer 



the channel estimation 

er is allocated for the data for 
about by better channel 



i 



j 20 

| 

bandwidth loss in contrast to the PSAM scheme which wil 
bandwidth, 



In general the BER is relatively srobust to the selection of £ and T|, and is 
consistent for other values of A* ( = AJ, A* ( = AJ, and diiferent values of Doppler 
spread. Optimum values selected are as B ff = 0.43 and Tf) 0 = 4/3 that minimize BER. 

With reference to PIG. 8, assuming that the SNR is known, the graph shows a 



plot of BER vs SNR for the prior art switched based blind 



scheme 801, and PSAM 802 



scheme, and the APSB scheme 803 and 804, using 6, - 0.43 and T\ Q = 4/3, at/ D r, = 
0.005. It is observed that when A„ = A L .= 15, plot 803, the SNR of the of the APSB 
10 scheme is 1 dB worse than the PSAM scheme. However, the BER performance of the 
APSB scheme is superior to the switched based Wind scheme. 

The present invention, as described, provides an added pilot semi-blind scheme 
that does not consume bandwidth, and is suitable for use jp a mobile communication 
system. 

IS This is accomplished by adding data and pilot data at the transmitter prior to 

transmission on a communication channel, and using a 2-JD Wiener filter to recover the 
transmission characteristics of the communication channe! using the pilot data. Then 
through an iterative process of estimating the communication channel and estimating 
the data, the communication channel can be estimated and the transmitted data 
20 recovered. 

The present invention provides a method and an applaratus for semi-blind 
communication channel estimation, which overcomes, or ; 
abovementioned problems of the prior art 

It will be appreciated that although only one particular embodiment of the 
25 invention has been described in detail, various modificatiens and improvements can be 
made by a person skilled in the art without departing from the scope of the present 
invention. 



always, have loss of 



;it least reduces the 



We claim; 



